1 HOK comparison of CASAL and Casal2 model configurations

This document compares the results of at least 2 CASAL model configurations (base and at least one sensitivity) and up to 6 Casal2 model configurations (3 BetaDiff and 3 ADOL-C).

The CASAL model sensitivity 1 has a smaller minimisation tolerance value than the CASAL base model (1e-9 vs. 1e-6).

The Casal2 ADOL-C and BetaDiff low tolerance models have a smaller tolerance value than the CASAL base model (1e-9 vs. 1e-6).

1.1 HOK model characteristics

The main characteristics of the Test Case HOK (hoki) CASAL model are:

  • one stock, ages 1 - 17
  • one area, although many characteristics are separated into “east” and “west”
  • years 1972 - 2018, projection years 2019 - 2023
  • five time steps: Oct_Nov, Dec_Mar, Apr_Jun, End_Jun, and Jul_Sep
  • five stock categories: west.sa, east.cr, west.cr, west.wc, and east.cs (labeled “stock.area”)
  • five migration processes, which have age-specific proportions specified
  • sex-specific natural mortality-at-age is double exponential
  • area-specific von Bertalanffy age-length relationship
  • area-specific length-weight relationship (\(W = aL^b\))
  • area-specific Beverton-Holt stock-recruitment relationships, with steepness (h) 0.75
  • ageing error is specified
  • surveys with double normal selectivity
  • six area-specific fisheries with double normal selectivity

Observation data include:

  • survey indices
  • survey proportions-at-age
  • fishery proportions-at-age

Parameters estimated include:

  • area-specific B0
  • parameters for the natural mortality-at-age ogive
  • catchability (q) for the surveys
  • parameters for the selectivity curves for the surveys
  • parameters for the selectivity curves for the fisheries
  • migration proportions by age
  • area-specific YCS (recruitment deviations)

1.2 R environment

## [1] "Wed Nov  3 11:47:17 2021"
## R version 3.6.0 (2019-04-26)
## Platform: x86_64-redhat-linux-gnu (64-bit)
## Running under: CentOS Linux 7 (Core)
## 
## Matrix products: default
## BLAS/LAPACK: /usr/lib64/R/lib/libRblas.so
## 
## locale:
##  [1] LC_CTYPE=en_NZ.UTF-8       LC_NUMERIC=C              
##  [3] LC_TIME=en_NZ.UTF-8        LC_COLLATE=en_NZ.UTF-8    
##  [5] LC_MONETARY=en_NZ.UTF-8    LC_MESSAGES=en_NZ.UTF-8   
##  [7] LC_PAPER=en_GB.UTF_8       LC_NAME=C                 
##  [9] LC_ADDRESS=C               LC_TELEPHONE=C            
## [11] LC_MEASUREMENT=en_NZ.UTF-8 LC_IDENTIFICATION=C       
## 
## attached base packages:
## [1] stats     graphics  grDevices utils     datasets  methods   base     
## 
## other attached packages:
##  [1] Casal2_21.10   casal_2.30     devtools_2.4.2 usethis_2.1.3  rlist_0.4.6.2 
##  [6] ggthemes_4.2.4 gridExtra_2.3  coda_0.19-4    ggmcmc_1.5.1.1 ggplot2_3.3.5 
## [11] tidyr_1.1.4    huxtable_5.4.0 dplyr_1.0.7    plyr_1.8.6    
## 
## loaded via a namespace (and not attached):
##  [1] Rcpp_1.0.7         lattice_0.20-45    prettyunits_1.1.1  ps_1.6.0          
##  [5] assertthat_0.2.1   rprojroot_2.0.2    digest_0.6.28      utf8_1.2.2        
##  [9] R6_2.5.1           evaluate_0.14      pillar_1.6.4       rlang_0.4.12      
## [13] data.table_1.14.2  callr_3.7.0        jquerylib_0.1.4    rmarkdown_2.11    
## [17] desc_1.4.0         stringr_1.4.0      munsell_0.5.0      compiler_3.6.0    
## [21] xfun_0.27          pkgconfig_2.0.3    pkgbuild_1.2.0     htmltools_0.5.2   
## [25] tidyselect_1.1.1   tibble_3.1.5       reshape_0.8.8      fansi_0.5.0       
## [29] crayon_1.4.2       withr_2.4.2        grid_3.6.0         jsonlite_1.7.2    
## [33] GGally_2.1.2       gtable_0.3.0       lifecycle_1.0.1    DBI_1.1.1         
## [37] magrittr_2.0.1     scales_1.1.1       cli_3.1.0          stringi_1.7.5     
## [41] cachem_1.0.6       fs_1.5.0           remotes_2.4.1      testthat_3.1.0    
## [45] bslib_0.3.1        ellipsis_0.3.2     generics_0.1.1     vctrs_0.3.8       
## [49] RColorBrewer_1.1-2 tools_3.6.0        glue_1.4.2         purrr_0.3.4       
## [53] processx_3.5.2     pkgload_1.2.3      fastmap_1.1.0      yaml_2.2.1        
## [57] colorspace_2.0-2   sessioninfo_1.2.0  memoise_2.0.0      knitr_1.36        
## [61] sass_0.4.0

1.3 CASAL and Casal2 model output

source('../../R-functions/report_read_in_CASAL_MPD_files.R')
source('../../R-functions/report_read_in_Casal2_MPD_files.R')

1.4 Tables

Tables of parameter estimates and objective function components for the CASAL and Casal2 model MPD results

CASAL parameter estimates
rownamesBase_ModelSensitivity_1Percent_Diff
q[CSacous].q               0.431200.431200.000
q[WCacous].q               0.300300.300200.033
q[CRsum].q               0.074310.074310.000
q[SAsum].q               0.059250.059250.000
q[SAaut].q               0.067170.067170.000
relative_abundance[CRsumbio].cv_process_error0.140500.140500.000
relative_abundance[SAsumbio].cv_process_error0.430600.430600.000
log_B0_total                 14.3400014.340000.000
B0_prop_stock1                 0.328700.328700.000
natural_mortality.all             0.296000.296000.000
selectivity[Wspsl].shift_a         -0.21740-0.217400.000
recruitment[E].YCS1            0.447800.447800.000
recruitment[E].YCS2            0.726800.726800.000
recruitment[E].YCS3            0.784500.784500.000
recruitment[E].YCS4            1.360001.360000.000
recruitment[E].YCS5            0.987700.98780-0.010
recruitment[E].YCS6            0.329300.329300.000
recruitment[E].YCS7            0.447200.447200.000
recruitment[E].YCS8            0.265600.265600.000
recruitment[E].YCS9            0.434200.434200.000
recruitment[E].YCS10           0.488000.488000.000
recruitment[E].YCS11           0.395700.395700.000
recruitment[E].YCS12           0.217900.217900.000
recruitment[E].YCS13           1.036001.036000.000
recruitment[E].YCS14           0.958900.958900.000
recruitment[E].YCS15           0.429000.429000.000
recruitment[E].YCS16           0.138100.138100.000
recruitment[E].YCS17           0.375000.375000.000
recruitment[E].YCS18           0.391700.391600.026
recruitment[E].YCS19           0.383100.383100.000
recruitment[E].YCS20           0.562900.562900.000
recruitment[E].YCS21           0.264000.264000.000
recruitment[E].YCS22           0.284500.284500.000
recruitment[E].YCS23           0.450000.450000.000
recruitment[E].YCS24           0.253100.253100.000
recruitment[E].YCS25           0.060000.060000.000
recruitment[E].YCS26           0.462200.462200.000
recruitment[E].YCS27           0.161600.161600.000
recruitment[E].YCS28           0.266900.266900.000
recruitment[E].YCS29           0.500400.500400.000
recruitment[E].YCS30           0.520000.520000.000
recruitment[E].YCS31           0.508800.508800.000
recruitment[E].YCS32           0.600600.600600.000
recruitment[E].YCS33           0.522100.522100.000
recruitment[E].YCS34           0.373000.373000.000
recruitment[E].YCS35           0.526800.526800.000
recruitment[E].YCS36           0.158600.158600.000
recruitment[E].YCS37           0.660600.660600.000
recruitment[E].YCS38           0.156000.156000.000
recruitment[E].YCS39           0.185800.185800.000
recruitment[E].YCS40           0.730600.730600.000
recruitment[E].YCS41           0.739100.739000.014
recruitment[E].YCS42           1.195001.195000.000
recruitment[W].YCS1            0.457000.457000.000
recruitment[W].YCS2            0.822200.822200.000
recruitment[W].YCS3            0.670000.670000.000
recruitment[W].YCS4            0.653900.653900.000
recruitment[W].YCS5            0.465900.465900.000
recruitment[W].YCS6            0.515100.515100.000
recruitment[W].YCS7            0.578200.578200.000
recruitment[W].YCS8            0.473300.473300.000
recruitment[W].YCS9            0.871900.871900.000
recruitment[W].YCS10           0.836400.83650-0.012
recruitment[W].YCS11           0.320200.320200.000
recruitment[W].YCS12           0.165000.165000.000
recruitment[W].YCS13           1.476001.476000.000
recruitment[W].YCS14           0.772700.77280-0.013
recruitment[W].YCS15           0.081350.081350.000
recruitment[W].YCS16           0.081580.081580.000
recruitment[W].YCS17           0.892800.892800.000
recruitment[W].YCS18           1.268001.268000.000
recruitment[W].YCS19           0.377200.377200.000
recruitment[W].YCS20           0.959400.959400.000
recruitment[W].YCS21           0.221100.221100.000
recruitment[W].YCS22           0.197200.197200.000
recruitment[W].YCS23           0.237900.23800-0.042
recruitment[W].YCS24           0.202900.202900.000
recruitment[W].YCS25           0.060000.060000.000
recruitment[W].YCS26           0.373000.373000.000
recruitment[W].YCS27           0.095480.095480.000
recruitment[W].YCS28           0.795600.795600.000
recruitment[W].YCS29           0.465200.465200.000
recruitment[W].YCS30           0.516600.516600.000
recruitment[W].YCS31           0.511300.511300.000
recruitment[W].YCS32           0.755400.755400.000
recruitment[W].YCS33           0.638900.638900.000
recruitment[W].YCS34           0.616900.616900.000
recruitment[W].YCS35           0.578400.578400.000
recruitment[W].YCS36           0.145500.145500.000
recruitment[W].YCS37           0.974800.974800.000
recruitment[W].YCS38           0.160200.160200.000
recruitment[W].YCS39           0.244400.244400.000
recruitment[W].YCS40           2.197002.197000.000
recruitment[W].YCS41           2.152002.152000.000
recruitment[W].YCS42           1.195001.195000.000
migration[Whome].rates_all1      0.327800.327800.000
migration[Whome].rates_all2      0.307800.307800.000
migration[Whome].rates_all3      0.476200.476200.000
migration[Whome].rates_all4      0.141800.141800.000
migration[Whome].rates_all5      0.468500.468400.021
migration[Whome].rates_all6      0.749800.75030-0.067
migration[Whome].rates_all7      0.794700.788400.793
migration[Whome].rates_all8      1.000001.000000.000
migration[Espmg].rates_all1      0.005660.005660.000
migration[Espmg].rates_all2      0.160200.160200.000
migration[Espmg].rates_all3      0.324100.324100.000
migration[Espmg].rates_all4      0.431200.431200.000
migration[Espmg].rates_all5      0.385500.385500.000
migration[Espmg].rates_all6      0.465400.465300.021
migration[Espmg].rates_all7      0.599000.599000.000
migration[Espmg].rates_all8      0.607800.607800.000
migration[Wspmg].rates_all1      0.106100.106100.000
migration[Wspmg].rates_all2      0.287700.287700.000
migration[Wspmg].rates_all3      0.270800.270800.000
migration[Wspmg].rates_all4      0.455300.455300.000
migration[Wspmg].rates_all5      0.484700.484700.000
migration[Wspmg].rates_all6      0.527100.527100.000
migration[Wspmg].rates_all7      0.660600.660600.000
migration[Wspmg].rates_all8      0.596300.596300.000
selectivity[Enspsl].all1            3.567003.567000.000
selectivity[Enspsl].all2            1.221001.221000.000
selectivity[Enspsl].all3            15.7100015.710000.000
selectivity[Wnspsl].all1            5.341005.341000.000
selectivity[Wnspsl].all2            2.839002.839000.000
selectivity[Wnspsl].all3            13.3800013.380000.000
selectivity[Espsl].all1            10.8800010.880000.000
selectivity[Espsl].all2            11.6800011.680000.000
selectivity[Espsl].all3            3.982003.982000.000
selectivity[CRsl].all1            1.000001.000000.000
selectivity[CRsl].all2            1.000001.000000.000
selectivity[CRsl].all3            9.300009.300000.000
selectivity[SAsl].all1            14.2500014.250000.000
selectivity[SAsl].all2            8.755008.75600-0.011
selectivity[SAsl].all3            1.315001.315000.000
Casal2 parameter estimates
rownamesbetadiff_casal_flags_onbetadiff_casal_flags_offbetadiff_casal_flags_on_low_toladolc_casal_flags_onadolc_casal_flags_offadolc_casal_flags_on_low_tol
catchability[CSacous].q               0.427860.427240.431160.427240.427240.42786
catchability[WCacous].q               0.312440.312840.300240.312840.312840.31244
catchability[CRsum].q               0.068110.068320.074310.068320.068320.06811
catchability[SAsum].q               0.067060.067170.059250.067170.067170.06706
catchability[SAaut].q               0.075630.075780.067170.075780.075780.07563
observation[CRsumbio].process_error   0.140890.140940.140480.140940.140940.14089
observation[SAsumbio].process_error   0.423810.423640.430610.423640.423640.42381
process[recruit_E].b0              547497.00000548395.00000555125.00000548396.00000548395.00000547497.00000
process[recruit_W].b0              1083360.000001082670.000001133680.000001082670.000001082670.000001083360.00000
process[Instant_mortality].m{west.sa}             0.312930.312710.295980.312710.312710.31293
time_varying[shifted_mu].a               -0.22959-0.22921-0.21744-0.22921-0.22921-0.22959
process[recruit_E].ycs_values{1975}0.445900.445380.447820.445380.445370.44589
process[recruit_E].ycs_values{1976}0.730520.729390.726800.729390.729390.73052
process[recruit_E].ycs_values{1977}0.793320.792310.784540.792310.792310.79332
process[recruit_E].ycs_values{1978}1.388271.386451.360011.386451.386451.38827
process[recruit_E].ycs_values{1979}1.003721.002650.987761.002651.002651.00372
process[recruit_E].ycs_values{1980}0.334620.334370.329320.334370.334370.33462
process[recruit_E].ycs_values{1981}0.455350.455020.447230.455020.455020.45535
process[recruit_E].ycs_values{1982}0.269890.269740.265560.269740.269740.26989
process[recruit_E].ycs_values{1983}0.441000.441060.434180.441060.441060.44100
process[recruit_E].ycs_values{1984}0.494230.494120.488040.494120.494120.49422
process[recruit_E].ycs_values{1985}0.400140.400100.395700.400100.400100.40014
process[recruit_E].ycs_values{1986}0.219820.219780.217910.219780.219780.21982
process[recruit_E].ycs_values{1987}1.040951.041311.036401.041311.041311.04095
process[recruit_E].ycs_values{1988}0.963150.963460.958910.963460.963460.96315
process[recruit_E].ycs_values{1989}0.431120.430950.428990.430950.430950.43112
process[recruit_E].ycs_values{1990}0.137980.137970.138100.137970.137970.13798
process[recruit_E].ycs_values{1991}0.373350.373850.375000.373850.373850.37335
process[recruit_E].ycs_values{1992}0.387570.388580.391640.388580.388580.38757
process[recruit_E].ycs_values{1993}0.379880.379950.383100.379950.379950.37988
process[recruit_E].ycs_values{1994}0.557260.557940.562850.557940.557940.55726
process[recruit_E].ycs_values{1995}0.261780.261780.263990.261780.261780.26178
process[recruit_E].ycs_values{1996}0.282480.282510.284470.282510.282510.28248
process[recruit_E].ycs_values{1997}0.447590.447810.450030.447810.447810.44759
process[recruit_E].ycs_values{1998}0.252490.252640.253070.252640.252640.25249
process[recruit_E].ycs_values{1999}0.060000.060000.060000.060000.060000.06000
process[recruit_E].ycs_values{2000}0.459870.460130.462200.460130.460130.45987
process[recruit_E].ycs_values{2001}0.160870.160860.161650.160860.160860.16087
process[recruit_E].ycs_values{2002}0.265600.266110.266890.266110.266110.26560
process[recruit_E].ycs_values{2003}0.499430.499590.500380.499590.499590.49943
process[recruit_E].ycs_values{2004}0.520310.520530.520030.520530.520530.52031
process[recruit_E].ycs_values{2005}0.507540.507640.508820.507640.507640.50754
process[recruit_E].ycs_values{2006}0.599860.600130.600550.600130.600130.59986
process[recruit_E].ycs_values{2007}0.520800.520870.522150.520870.520870.52080
process[recruit_E].ycs_values{2008}0.370460.371150.373000.371150.371150.37047
process[recruit_E].ycs_values{2009}0.527860.527900.526810.527890.527900.52786
process[recruit_E].ycs_values{2010}0.156900.156950.158640.156950.156950.15690
process[recruit_E].ycs_values{2011}0.660290.660470.660580.660470.660470.66029
process[recruit_E].ycs_values{2012}0.155580.155540.156000.155540.155540.15558
process[recruit_E].ycs_values{2013}0.185180.185180.185840.185180.185180.18518
process[recruit_E].ycs_values{2014}0.728080.728440.730620.728440.728440.72808
process[recruit_E].ycs_values{2015}0.752300.755860.739000.755870.755860.75231
process[recruit_E].ycs_values{2016}1.189391.188861.195151.188861.188861.18939
process[recruit_W].ycs_values{1975}0.458720.458690.457020.458690.458690.45872
process[recruit_W].ycs_values{1976}0.814330.814220.822180.814220.814220.81433
process[recruit_W].ycs_values{1977}0.684940.685150.670000.685140.685150.68494
process[recruit_W].ycs_values{1978}0.669050.669370.653900.669380.669370.66905
process[recruit_W].ycs_values{1979}0.472760.473070.465900.473070.473070.47276
process[recruit_W].ycs_values{1980}0.520520.520640.515120.520640.520640.52052
process[recruit_W].ycs_values{1981}0.583280.583440.578240.583440.583440.58328
process[recruit_W].ycs_values{1982}0.474400.474620.473280.474620.474620.47440
process[recruit_W].ycs_values{1983}0.873750.873090.871920.873090.873090.87375
process[recruit_W].ycs_values{1984}0.837380.836330.836460.836330.836330.83738
process[recruit_W].ycs_values{1985}0.319750.319420.320240.319420.319420.31975
process[recruit_W].ycs_values{1986}0.164840.164930.164970.164930.164930.16484
process[recruit_W].ycs_values{1987}1.480661.480071.475551.480071.480071.48066
process[recruit_W].ycs_values{1988}0.775300.774540.772760.774540.774540.77530
process[recruit_W].ycs_values{1989}0.081470.081480.081340.081480.081480.08147
process[recruit_W].ycs_values{1990}0.081800.081830.081580.081830.081830.08180
process[recruit_W].ycs_values{1991}0.891890.891970.892800.891970.891970.89189
process[recruit_W].ycs_values{1992}1.267831.267571.267781.267571.267571.26783
process[recruit_W].ycs_values{1993}0.376480.376520.377220.376520.376520.37648
process[recruit_W].ycs_values{1994}0.955570.955790.959390.955790.955790.95557
process[recruit_W].ycs_values{1995}0.220150.220260.221130.220260.220260.22015
process[recruit_W].ycs_values{1996}0.196260.196320.197180.196320.196320.19626
process[recruit_W].ycs_values{1997}0.236780.236430.237950.236430.236430.23678
process[recruit_W].ycs_values{1998}0.200610.200410.202920.200410.200410.20061
process[recruit_W].ycs_values{1999}0.060000.060000.060000.060000.060000.06000
process[recruit_W].ycs_values{2000}0.371330.370890.373040.370890.370890.37133
process[recruit_W].ycs_values{2001}0.095100.095060.095480.095060.095060.09510
process[recruit_W].ycs_values{2002}0.795760.795530.795580.795530.795530.79576
process[recruit_W].ycs_values{2003}0.465880.465700.465170.465700.465700.46588
process[recruit_W].ycs_values{2004}0.517840.517610.516620.517610.517610.51784
process[recruit_W].ycs_values{2005}0.513140.513050.511260.513050.513050.51314
process[recruit_W].ycs_values{2006}0.754870.754750.755440.754750.754750.75487
process[recruit_W].ycs_values{2007}0.636760.636920.638910.636920.636920.63676
process[recruit_W].ycs_values{2008}0.616140.615640.616860.615640.615640.61614
process[recruit_W].ycs_values{2009}0.572920.573050.578410.573050.573050.57292
process[recruit_W].ycs_values{2010}0.146260.146180.145470.146180.146180.14626
process[recruit_W].ycs_values{2011}0.966690.967080.974780.967080.967080.96669
process[recruit_W].ycs_values{2012}0.158970.159090.160160.159090.159090.15897
process[recruit_W].ycs_values{2013}0.243890.244020.244440.244020.244020.24389
process[recruit_W].ycs_values{2014}2.184742.186512.196572.186512.186512.18474
process[recruit_W].ycs_values{2015}2.149692.147582.152042.147572.147582.14969
process[recruit_W].ycs_values{2016}1.189621.189091.195391.189091.189091.18962
selectivity[sel_Whome].v{1}            0.321360.321930.327830.321930.321930.32136
selectivity[sel_Whome].v{2}            0.283860.284910.307780.284910.284910.28386
selectivity[sel_Whome].v{3}            0.457750.459720.476200.459730.459720.45775
selectivity[sel_Whome].v{4}            0.135020.135490.141840.135490.135490.13502
selectivity[sel_Whome].v{5}            0.431690.437820.468360.437830.437820.43169
selectivity[sel_Whome].v{6}            0.718920.708210.750260.708210.708200.71891
selectivity[sel_Whome].v{7}            0.759540.999870.788410.999900.999880.75960
selectivity[sel_Espmg].v{1}            0.005640.005640.005660.005640.005640.00564
selectivity[sel_Espmg].v{2}            0.159010.159120.160210.159120.159120.15901
selectivity[sel_Espmg].v{3}            0.320560.320820.324110.320820.320820.32056
selectivity[sel_Espmg].v{4}            0.425130.425470.431230.425470.425470.42513
selectivity[sel_Espmg].v{5}            0.379270.379530.385460.379530.379530.37927
selectivity[sel_Espmg].v{6}            0.458470.458760.465350.458760.458760.45847
selectivity[sel_Espmg].v{7}            0.590440.590810.599020.590810.590810.59044
selectivity[sel_Espmg].v{8}            0.610640.611000.607780.611000.611000.61064
selectivity[sel_Wspmg].v{1}            0.105530.105560.106140.105560.105560.10553
selectivity[sel_Wspmg].v{2}            0.291390.291360.287740.291360.291360.29139
selectivity[sel_Wspmg].v{3}            0.270660.270680.270760.270680.270680.27066
selectivity[sel_Wspmg].v{4}            0.455200.455270.455250.455270.455270.45520
selectivity[sel_Wspmg].v{5}            0.484340.484210.484660.484210.484210.48434
selectivity[sel_Wspmg].v{6}            0.523980.524960.527110.524960.524960.52398
selectivity[sel_Wspmg].v{7}            0.655910.652980.660620.652980.652980.65591
selectivity[sel_Wspmg].v{8}            0.601010.601530.596270.601530.601530.60101
selectivity[Enspsl].mu              3.572273.573183.567163.573183.573183.57227
selectivity[Enspsl].sigma_l         1.222061.222451.221081.222451.222451.22206
selectivity[Enspsl].sigma_r         49.9989050.0000015.7137050.0000049.9999049.99990
selectivity[Wnspsl].mu              5.329515.328405.341175.328405.328405.32952
selectivity[Wnspsl].sigma_l         2.778452.778152.839342.778152.778152.77846
selectivity[Wnspsl].sigma_r         49.9949049.9999013.3789049.9998050.0000049.99930
selectivity[Espsl].mu              11.2546011.2567010.8788011.2567011.2567011.25470
selectivity[Espsl].sigma_l         10.8594010.8674011.6775010.8674010.8674010.85940
selectivity[Espsl].sigma_r         4.029124.026253.981654.026254.026254.02912
selectivity[CRsl].mu              1.000001.000001.000001.000001.000001.00000
selectivity[CRsl].sigma_l         46.5212046.5222046.0231046.5222046.5219046.52180
selectivity[CRsl].sigma_r         10.0878010.092809.3004210.0928010.0928010.08780
selectivity[SAsl].mu              14.8648014.8685014.2523014.8685014.8685014.86480
selectivity[SAsl].sigma_l         8.692818.694688.755588.694698.694688.69281
selectivity[SAsl].sigma_r         1.000001.000001.315191.000001.000001.00000
Casal2 parameter estimates: Percent Difference from betadiff_casal_flags_on
rownamesbetadiff_casal_flags_onbetadiff_casal_flags_offbetadiff_casal_flags_on_low_toladolc_casal_flags_onadolc_casal_flags_offadolc_casal_flags_on_low_tol
catchability[CSacous].q               0.0000.145-0.7720.1450.1450.000
catchability[WCacous].q               0.000-0.1283.904-0.128-0.1280.000
catchability[CRsum].q               0.000-0.299-9.090-0.299-0.299-0.000
catchability[SAsum].q               0.000-0.17211.645-0.172-0.172-0.000
catchability[SAaut].q               0.000-0.18911.196-0.189-0.189-0.000
observation[CRsumbio].process_error   0.000-0.0350.292-0.035-0.0350.000
observation[SAsumbio].process_error   0.0000.041-1.6040.0410.0400.000
process[recruit_E].b0              0.000-0.164-1.393-0.164-0.1640.000
process[recruit_W].b0              0.0000.064-4.6450.0640.0640.000
process[Instant_mortality].m{west.sa}             0.0000.0705.4150.0700.0700.000
time_varying[shifted_mu].a               0.0000.1665.2930.1660.1660.000
process[recruit_E].ycs_values{1975}0.0000.117-0.4310.1170.1170.000
process[recruit_E].ycs_values{1976}0.0000.1550.5100.1550.155-0.000
process[recruit_E].ycs_values{1977}0.0000.1281.1080.1280.1280.000
process[recruit_E].ycs_values{1978}0.0000.1312.0360.1310.1310.000
process[recruit_E].ycs_values{1979}0.0000.1071.5900.1070.1070.000
process[recruit_E].ycs_values{1980}0.0000.0741.5840.0740.0740.000
process[recruit_E].ycs_values{1981}0.0000.0721.7830.0720.0720.000
process[recruit_E].ycs_values{1982}0.0000.0551.6050.0550.0550.000
process[recruit_E].ycs_values{1983}0.000-0.0121.547-0.012-0.0120.000
process[recruit_E].ycs_values{1984}0.0000.0221.2510.0220.0220.000
process[recruit_E].ycs_values{1985}0.0000.0091.1110.0090.009-0.000
process[recruit_E].ycs_values{1986}0.0000.0220.8720.0220.0220.000
process[recruit_E].ycs_values{1987}0.000-0.0350.437-0.035-0.0350.000
process[recruit_E].ycs_values{1988}0.000-0.0320.440-0.032-0.0320.000
process[recruit_E].ycs_values{1989}0.0000.0400.4950.0400.040-0.000
process[recruit_E].ycs_values{1990}0.0000.009-0.0910.0090.0090.000
process[recruit_E].ycs_values{1991}0.000-0.136-0.443-0.135-0.1360.000
process[recruit_E].ycs_values{1992}0.000-0.260-1.050-0.260-0.2600.000
process[recruit_E].ycs_values{1993}0.000-0.016-0.847-0.016-0.0160.000
process[recruit_E].ycs_values{1994}0.000-0.123-1.003-0.123-0.123-0.000
process[recruit_E].ycs_values{1995}0.0000.002-0.8420.0020.0020.000
process[recruit_E].ycs_values{1996}0.000-0.008-0.704-0.008-0.0080.000
process[recruit_E].ycs_values{1997}0.000-0.049-0.546-0.049-0.049-0.000
process[recruit_E].ycs_values{1998}0.000-0.062-0.232-0.062-0.062-0.000
process[recruit_E].ycs_values{1999}0.0000.0000.0000.0000.0000.000
process[recruit_E].ycs_values{2000}0.000-0.057-0.506-0.057-0.0570.000
process[recruit_E].ycs_values{2001}0.0000.007-0.4820.0070.0070.000
process[recruit_E].ycs_values{2002}0.000-0.191-0.486-0.191-0.1910.000
process[recruit_E].ycs_values{2003}0.000-0.031-0.190-0.031-0.0310.000
process[recruit_E].ycs_values{2004}0.000-0.0430.054-0.043-0.0430.000
process[recruit_E].ycs_values{2005}0.000-0.021-0.254-0.021-0.0210.000
process[recruit_E].ycs_values{2006}0.000-0.046-0.116-0.046-0.0460.000
process[recruit_E].ycs_values{2007}0.000-0.014-0.259-0.014-0.0140.000
process[recruit_E].ycs_values{2008}0.000-0.184-0.686-0.184-0.184-0.000
process[recruit_E].ycs_values{2009}0.000-0.0070.199-0.007-0.0070.000
process[recruit_E].ycs_values{2010}0.000-0.033-1.107-0.033-0.0330.000
process[recruit_E].ycs_values{2011}0.000-0.028-0.044-0.028-0.0280.000
process[recruit_E].ycs_values{2012}0.0000.030-0.2640.0300.0300.000
process[recruit_E].ycs_values{2013}0.000-0.002-0.357-0.002-0.0020.000
process[recruit_E].ycs_values{2014}0.000-0.049-0.349-0.049-0.0490.000
process[recruit_E].ycs_values{2015}0.000-0.4731.768-0.474-0.473-0.000
process[recruit_E].ycs_values{2016}0.0000.045-0.4840.0450.0450.000
process[recruit_W].ycs_values{1975}0.0000.0050.3700.0050.006-0.000
process[recruit_W].ycs_values{1976}0.0000.014-0.9630.0140.0140.000
process[recruit_W].ycs_values{1977}0.000-0.0302.181-0.030-0.0300.000
process[recruit_W].ycs_values{1978}0.000-0.0482.265-0.049-0.048-0.001
process[recruit_W].ycs_values{1979}0.000-0.0671.451-0.067-0.0670.000
process[recruit_W].ycs_values{1980}0.000-0.0231.038-0.023-0.0230.000
process[recruit_W].ycs_values{1981}0.000-0.0270.865-0.027-0.027-0.000
process[recruit_W].ycs_values{1982}0.000-0.0470.237-0.048-0.047-0.000
process[recruit_W].ycs_values{1983}0.0000.0760.2090.0760.0760.000
process[recruit_W].ycs_values{1984}0.0000.1250.1090.1250.1250.000
process[recruit_W].ycs_values{1985}0.0000.103-0.1530.1030.1030.000
process[recruit_W].ycs_values{1986}0.000-0.054-0.075-0.054-0.0540.001
process[recruit_W].ycs_values{1987}0.0000.0400.3450.0400.0400.000
process[recruit_W].ycs_values{1988}0.0000.0980.3270.0990.0990.000
process[recruit_W].ycs_values{1989}0.000-0.0030.158-0.003-0.0030.000
process[recruit_W].ycs_values{1990}0.000-0.0370.267-0.037-0.0370.000
process[recruit_W].ycs_values{1991}0.000-0.009-0.102-0.009-0.009-0.000
process[recruit_W].ycs_values{1992}0.0000.0210.0040.0210.0210.000
process[recruit_W].ycs_values{1993}0.000-0.011-0.196-0.011-0.0110.000
process[recruit_W].ycs_values{1994}0.000-0.023-0.400-0.023-0.0230.000
process[recruit_W].ycs_values{1995}0.000-0.054-0.446-0.054-0.0540.000
process[recruit_W].ycs_values{1996}0.000-0.028-0.468-0.028-0.0280.000
process[recruit_W].ycs_values{1997}0.0000.147-0.4950.1470.1470.000
process[recruit_W].ycs_values{1998}0.0000.100-1.1500.1000.1000.000
process[recruit_W].ycs_values{1999}0.0000.0000.0000.0000.0000.000
process[recruit_W].ycs_values{2000}0.0000.117-0.4620.1180.1180.000
process[recruit_W].ycs_values{2001}0.0000.038-0.4090.0380.0380.000
process[recruit_W].ycs_values{2002}0.0000.0290.0220.0290.0290.000
process[recruit_W].ycs_values{2003}0.0000.0390.1530.0390.0390.000
process[recruit_W].ycs_values{2004}0.0000.0450.2350.0450.0450.000
process[recruit_W].ycs_values{2005}0.0000.0170.3660.0170.0170.000
process[recruit_W].ycs_values{2006}0.0000.015-0.0760.0150.0150.000
process[recruit_W].ycs_values{2007}0.000-0.025-0.338-0.025-0.0250.000
process[recruit_W].ycs_values{2008}0.0000.081-0.1180.0810.0810.000
process[recruit_W].ycs_values{2009}0.000-0.024-0.958-0.024-0.0240.000
process[recruit_W].ycs_values{2010}0.0000.0550.5380.0550.0550.000
process[recruit_W].ycs_values{2011}0.000-0.041-0.837-0.041-0.0410.000
process[recruit_W].ycs_values{2012}0.000-0.075-0.750-0.075-0.0750.000
process[recruit_W].ycs_values{2013}0.000-0.052-0.224-0.052-0.0520.000
process[recruit_W].ycs_values{2014}0.000-0.081-0.541-0.081-0.0810.000
process[recruit_W].ycs_values{2015}0.0000.098-0.1090.0990.0980.000
process[recruit_W].ycs_values{2016}0.0000.045-0.4850.0450.0450.000
selectivity[sel_Whome].v{1}            0.000-0.177-2.014-0.177-0.1770.000
selectivity[sel_Whome].v{2}            0.000-0.371-8.428-0.371-0.3710.000
selectivity[sel_Whome].v{3}            0.000-0.431-4.031-0.432-0.4310.000
selectivity[sel_Whome].v{4}            0.000-0.344-5.047-0.342-0.3430.001
selectivity[sel_Whome].v{5}            0.000-1.421-8.495-1.422-1.421-0.000
selectivity[sel_Whome].v{6}            0.0001.489-4.3601.4891.4900.001
selectivity[sel_Whome].v{7}            0.000-31.641-3.801-31.645-31.642-0.008
selectivity[sel_Espmg].v{1}            0.000-0.081-0.412-0.081-0.0810.000
selectivity[sel_Espmg].v{2}            0.000-0.073-0.757-0.073-0.0730.000
selectivity[sel_Espmg].v{3}            0.000-0.078-1.107-0.078-0.0780.000
selectivity[sel_Espmg].v{4}            0.000-0.079-1.434-0.079-0.0790.000
selectivity[sel_Espmg].v{5}            0.000-0.069-1.633-0.069-0.0690.000
selectivity[sel_Espmg].v{6}            0.000-0.063-1.500-0.063-0.0630.000
selectivity[sel_Espmg].v{7}            0.000-0.062-1.452-0.062-0.0620.000
selectivity[sel_Espmg].v{8}            0.000-0.0600.468-0.060-0.0600.000
selectivity[sel_Wspmg].v{1}            0.000-0.033-0.586-0.033-0.0330.000
selectivity[sel_Wspmg].v{2}            0.0000.0101.2520.0100.0100.000
selectivity[sel_Wspmg].v{3}            0.000-0.008-0.037-0.008-0.0080.000
selectivity[sel_Wspmg].v{4}            0.000-0.016-0.012-0.017-0.0170.000
selectivity[sel_Wspmg].v{5}            0.0000.026-0.0660.0260.0260.000
selectivity[sel_Wspmg].v{6}            0.000-0.187-0.597-0.187-0.1870.000
selectivity[sel_Wspmg].v{7}            0.0000.447-0.7180.4470.4470.000
selectivity[sel_Wspmg].v{8}            0.000-0.0880.788-0.088-0.0880.000
selectivity[Enspsl].mu              0.000-0.0250.143-0.025-0.0250.000
selectivity[Enspsl].sigma_l         0.000-0.0320.080-0.032-0.0320.000
selectivity[Enspsl].sigma_r         0.000-0.00268.572-0.002-0.002-0.002
selectivity[Wnspsl].mu              0.0000.021-0.2190.0210.021-0.000
selectivity[Wnspsl].sigma_l         0.0000.011-2.1920.0110.011-0.000
selectivity[Wnspsl].sigma_r         0.000-0.01073.239-0.010-0.010-0.009
selectivity[Espsl].mu              0.000-0.0193.339-0.019-0.019-0.001
selectivity[Espsl].sigma_l         0.000-0.074-7.534-0.074-0.0740.000
selectivity[Espsl].sigma_r         0.0000.0711.1780.0710.0710.000
selectivity[CRsl].mu              0.0000.0000.0000.0000.0000.000
selectivity[CRsl].sigma_l         0.000-0.0021.071-0.002-0.002-0.001
selectivity[CRsl].sigma_r         0.000-0.0507.805-0.050-0.0500.000
selectivity[SAsl].mu              0.000-0.0254.120-0.025-0.0250.000
selectivity[SAsl].sigma_l         0.000-0.022-0.722-0.022-0.0220.000
selectivity[SAsl].sigma_r         0.0000.000-31.5190.0000.0000.000

CASAL objective function component values
ComponentBase_ModelSensitivity_1
CRsumbio-31.030-31.030
CSacous-9.835-9.836
SAautbio-3.984-3.984
SAsumbio-6.146-6.146
WCacous-5.916-5.916
CRsumage323.800323.800
SAautage35.12035.120
SAsumage221.100221.100
EnspOLF24.28024.280
Enspage256.100256.100
Espage581.500581.500
WnspOLF70.25070.250
Wnspage157.000157.000
Wspage364.800364.800
pspawn-13.130-13.130
pspawn_1993-6.015-6.015
prior_on_q[CSacous].q-0.839-0.839
prior_on_q[WCacous].q-1.202-1.202
prior_on_q[CRsum].q-2.207-2.207
prior_on_q[SAsum].q-1.557-1.557
prior_on_q[SAaut].q-1.762-1.762
prior_on_relative_abundance[CRsumbio].cv_process_error0.0000.000
prior_on_relative_abundance[SAsumbio].cv_process_error0.0000.000
prior_on_log_B0_total0.0000.000
prior_on_B0_prop_stock112.56012.560
prior_on_natural_mortality.all-1.217-1.217
prior_on_selectivity[Wspsl].shift_a0.3780.378
prior_on_recruitment[E].YCS-13.750-13.750
prior_on_recruitment[W].YCS-1.785-1.785
prior_on_migration[Whome].rates_all0.0000.000
prior_on_migration[Espmg].rates_all0.0000.000
prior_on_migration[Wspmg].rates_all0.0000.000
prior_on_selectivity[Enspsl].all0.0000.000
prior_on_selectivity[Wnspsl].all0.0000.000
prior_on_selectivity[Espsl].all0.0000.000
prior_on_selectivity[CRsl].all0.0000.000
prior_on_selectivity[SAsl].all0.0000.000
Umax_Ensp10.0000.000
Umax_Ensp20.0000.000
Umax_Wnsp10.0000.000
Umax_Wnsp20.0000.000
Umax_Esp0.0000.000
Umax_Wsp0.0000.000
YCS.eq.160.0000.000
sp.migr0.0130.013
Total1947.0001947.000
Casal2 objective function component values
rownamesbetadiff_casal_flags_onbetadiff_casal_flags_offbetadiff_casal_flags_on_low_toladolc_casal_flags_onadolc_casal_flags_offadolc_casal_flags_on_low_tol
observation->CSacousbio-1991                 2.3442.3422.4042.3422.3422.344
observation->CSacousbio-1993                 -0.577-0.577-0.588-0.577-0.577-0.577
observation->CSacousbio-1994                 -0.016-0.016-0.030-0.016-0.016-0.016
observation->CSacousbio-1995                 -0.562-0.562-0.567-0.562-0.562-0.562
observation->CSacousbio-1996                 -0.027-0.0270.013-0.026-0.027-0.027
observation->CSacousbio-1997                 -0.912-0.912-0.899-0.912-0.912-0.912
observation->CSacousbio-1998                 -0.055-0.054-0.004-0.054-0.054-0.055
observation->CSacousbio-1999                 -1.025-1.025-1.012-1.025-1.025-1.025
observation->CSacousbio-2001                 -1.217-1.217-1.211-1.217-1.217-1.217
observation->CSacousbio-2002                 -0.297-0.297-0.338-0.297-0.297-0.297
observation->CSacousbio-2003                 -0.845-0.845-0.863-0.845-0.845-0.845
observation->CSacousbio-2005                 -0.929-0.929-0.927-0.929-0.929-0.929
observation->CSacousbio-2006                 -0.939-0.939-0.950-0.939-0.939-0.939
observation->CSacousbio-2007                 -0.623-0.623-0.600-0.623-0.623-0.623
observation->CSacousbio-2008                 -1.094-1.094-1.125-1.094-1.094-1.094
observation->CSacousbio-2009                 -0.186-0.185-0.110-0.185-0.185-0.186
observation->CSacousbio-2011                 -1.021-1.021-1.000-1.021-1.021-1.021
observation->CSacousbio-2013                 -0.892-0.892-0.855-0.892-0.892-0.892
observation->CSacousbio-2015                 -0.191-0.190-0.150-0.190-0.190-0.191
observation->CSacousbio-2017                 -1.020-1.020-1.021-1.020-1.020-1.020
observation->WCacousbio-1988                 -0.588-0.588-0.588-0.588-0.588-0.588
observation->WCacousbio-1989                 -0.184-0.183-0.168-0.183-0.183-0.184
observation->WCacousbio-1990                 -0.331-0.332-0.327-0.332-0.332-0.331
observation->WCacousbio-1991                 -0.425-0.425-0.425-0.425-0.425-0.425
observation->WCacousbio-1992                 -0.744-0.744-0.744-0.744-0.744-0.744
observation->WCacousbio-1993                 -0.371-0.369-0.374-0.369-0.369-0.371
observation->WCacousbio-1997                 -0.055-0.054-0.043-0.054-0.054-0.055
observation->WCacousbio-2000                 -1.264-1.265-1.272-1.265-1.265-1.264
observation->WCacousbio-2012                 -0.906-0.906-0.896-0.906-0.906-0.906
observation->WCacousbio-2013                 -1.079-1.079-1.079-1.079-1.079-1.079
observation->CRsumbio-1992                 -1.634-1.634-1.647-1.634-1.634-1.634
observation->CRsumbio-1993                 0.0410.0400.0430.0400.0400.041
observation->CRsumbio-1994                 -1.761-1.760-1.760-1.760-1.760-1.761
observation->CRsumbio-1995                 -1.431-1.425-1.406-1.425-1.425-1.431
observation->CRsumbio-1996                 -1.691-1.690-1.699-1.690-1.690-1.691
observation->CRsumbio-1997                 -0.609-0.607-0.635-0.607-0.607-0.609
observation->CRsumbio-1998                 -1.205-1.203-1.210-1.203-1.203-1.205
observation->CRsumbio-1999                 -0.950-0.953-0.919-0.953-0.953-0.950
observation->CRsumbio-2000                 -1.629-1.629-1.632-1.629-1.629-1.629
observation->CRsumbio-2001                 -1.773-1.773-1.775-1.773-1.773-1.773
observation->CRsumbio-2002                 -1.205-1.209-1.206-1.209-1.209-1.205
observation->CRsumbio-2003                 -1.748-1.749-1.744-1.749-1.749-1.748
observation->CRsumbio-2004                 -0.821-0.818-0.803-0.818-0.818-0.821
observation->CRsumbio-2005                 -1.470-1.471-1.477-1.471-1.471-1.470
observation->CRsumbio-2006                 -1.314-1.314-1.306-1.314-1.314-1.314
observation->CRsumbio-2007                 -0.726-0.727-0.756-0.727-0.727-0.726
observation->CRsumbio-2008                 -0.376-0.373-0.402-0.373-0.373-0.376
observation->CRsumbio-2009                 -0.587-0.586-0.559-0.586-0.586-0.587
observation->CRsumbio-2010                 -1.350-1.349-1.359-1.349-1.349-1.350
observation->CRsumbio-2011                 -1.094-1.096-1.105-1.096-1.096-1.094
observation->CRsumbio-2012                 -1.438-1.439-1.448-1.439-1.439-1.438
observation->CRsumbio-2013                 -1.494-1.494-1.491-1.494-1.494-1.494
observation->CRsumbio-2014                 -1.768-1.768-1.770-1.768-1.768-1.768
observation->CRsumbio-2016                 -1.411-1.410-1.411-1.410-1.410-1.411
observation->CRsumbio-2018                 -1.552-1.552-1.554-1.552-1.552-1.552
observation->SAsumbio-1992                 -0.596-0.597-0.588-0.597-0.597-0.596
observation->SAsumbio-1993                 -0.314-0.314-0.309-0.314-0.314-0.314
observation->SAsumbio-1994                 0.1940.1960.1860.1960.1960.194
observation->SAsumbio-2001                 -0.834-0.835-0.822-0.835-0.835-0.834
observation->SAsumbio-2002                 -0.794-0.794-0.777-0.794-0.794-0.794
observation->SAsumbio-2003                 -0.836-0.837-0.825-0.837-0.837-0.836
observation->SAsumbio-2004                 0.7350.7320.7340.7320.7320.735
observation->SAsumbio-2005                 -0.255-0.258-0.218-0.258-0.258-0.255
observation->SAsumbio-2006                 -0.598-0.600-0.558-0.600-0.600-0.598
observation->SAsumbio-2007                 0.6370.6360.6540.6360.6360.637
observation->SAsumbio-2008                 -0.463-0.461-0.462-0.461-0.461-0.463
observation->SAsumbio-2009                 -0.584-0.585-0.583-0.585-0.585-0.584
observation->SAsumbio-2010                 -0.146-0.143-0.147-0.143-0.143-0.146
observation->SAsumbio-2012                 -0.848-0.848-0.835-0.848-0.848-0.848
observation->SAsumbio-2013                 -0.796-0.797-0.786-0.797-0.797-0.796
observation->SAsumbio-2015                 -0.201-0.200-0.189-0.200-0.200-0.201
observation->SAsumbio-2017                 -0.638-0.638-0.620-0.638-0.638-0.638
observation->SAautbio-1992                 -1.460-1.459-1.450-1.459-1.459-1.460
observation->SAautbio-1996                 -1.232-1.233-1.250-1.233-1.233-1.232
observation->SAautbio-1998                 -1.306-1.308-1.284-1.308-1.308-1.306
observation->Espage-1988                 19.79019.79319.64019.79319.79319.790
observation->Espage-1989                 15.85715.85815.86215.85815.85815.857
observation->Espage-1990                 21.48621.48821.44621.48821.48821.486
observation->Espage-1991                 20.64020.64120.61220.64120.64120.640
observation->Espage-1992                 24.09224.08924.02024.08924.08924.092
observation->Espage-1993                 19.30419.30519.32819.30519.30519.304
observation->Espage-1994                 20.52220.52820.50320.52820.52820.522
observation->Espage-1995                 22.53022.54122.55922.54122.54122.530
observation->Espage-1996                 18.59018.59018.57718.59018.59018.590
observation->Espage-1997                 20.97820.97420.96420.97420.97420.978
observation->Espage-1998                 20.57920.58720.61120.58720.58720.579
observation->Espage-1999                 22.87522.88622.86222.88622.88622.875
observation->Espage-2000                 23.32123.32123.31423.32123.32123.321
observation->Espage-2001                 20.73520.73520.77320.73520.73520.735
observation->Espage-2002                 26.24326.24826.28726.24826.24826.243
observation->Espage-2003                 19.11819.12019.12019.12019.12019.118
observation->Espage-2004                 22.69822.69222.64222.69222.69222.698
observation->Espage-2005                 26.12026.12226.07026.12226.12226.120
observation->Espage-2006                 23.24523.24223.14923.24223.24223.245
observation->Espage-2007                 25.27725.27425.21025.27425.27425.277
observation->Espage-2008                 22.16622.16722.12922.16722.16722.166
observation->Espage-2009                 18.69118.69118.69018.69118.69118.691
observation->Espage-2010                 21.85621.85521.85121.85521.85521.856
observation->Espage-2014                 20.61420.61620.64820.61620.61620.614
observation->Espage-2015                 21.82021.82321.87021.82321.82321.820
observation->Espage-2016                 21.20321.20421.20421.20421.20421.203
observation->Espage-2017                 21.57821.57521.56321.57521.57521.578
observation->Wspage-1988                 15.72915.72415.75315.72415.72415.729
observation->Wspage-1989                 11.55811.55511.56611.55511.55511.558
observation->Wspage-1990                 11.75011.75211.79011.75211.75211.750
observation->Wspage-1991                 13.03713.03613.05113.03613.03613.037
observation->Wspage-1992                 10.17910.17810.17010.17810.17810.179
observation->Wspage-1993                 10.57110.56910.57910.56910.56910.571
observation->Wspage-1994                 9.9189.9189.9249.9189.9189.918
observation->Wspage-1995                 13.39913.40313.39613.40313.40313.399
observation->Wspage-1996                 11.84111.84011.84411.84011.84011.841
observation->Wspage-1997                 11.23311.23011.24311.23011.23011.233
observation->Wspage-1998                 11.76111.76111.76211.76111.76111.761
observation->Wspage-1999                 11.97611.97511.98111.97511.97511.976
observation->Wspage-2000                 11.25611.25611.26011.25611.25611.256
observation->Wspage-2001                 12.18112.17912.20112.17912.17912.181
observation->Wspage-2002                 9.7329.7319.7459.7319.7319.732
observation->Wspage-2003                 9.5559.5559.5569.5559.5559.555
observation->Wspage-2004                 14.34214.34014.30414.34014.34014.342
observation->Wspage-2005                 14.93414.93214.97714.93214.93214.934
observation->Wspage-2006                 13.81713.81713.87213.81713.81713.817
observation->Wspage-2007                 11.50011.49811.51611.49811.49811.500
observation->Wspage-2008                 12.19612.19512.19412.19512.19512.196
observation->Wspage-2009                 12.63512.63112.59412.63112.63112.635
observation->Wspage-2010                 11.31711.31611.29911.31611.31611.317
observation->Wspage-2011                 11.15511.15411.13711.15411.15411.155
observation->Wspage-2012                 11.86611.86711.86011.86711.86711.866
observation->Wspage-2013                 11.04311.04111.06511.04111.04111.043
observation->Wspage-2014                 12.94412.94412.94112.94412.94412.944
observation->Wspage-2015                 13.13613.13413.10113.13413.13413.136
observation->Wspage-2016                 14.26814.27214.27314.27214.27214.268
observation->Wspage-2017                 13.74613.74513.79713.74513.74513.746
observation->EnspOLF-1992                 4.3964.3974.3884.3974.3974.396
observation->EnspOLF-1994                 6.6456.6476.6786.6476.6476.645
observation->EnspOLF-1996                 7.8997.8977.9067.8977.8977.899
observation->EnspOLF-1998                 5.3015.3015.3035.3015.3015.301
observation->Enspage-1999                 14.33114.33214.28714.33214.33214.331
observation->Enspage-2000                 13.97313.97513.96013.97513.97513.973
observation->Enspage-2001                 13.95913.97913.92913.97913.97913.959
observation->Enspage-2002                 12.08712.08112.07812.08112.08112.087
observation->Enspage-2003                 9.7029.7009.7069.7009.7009.702
observation->Enspage-2004                 12.97212.96612.98612.96612.96612.972
observation->Enspage-2005                 11.00911.01211.00911.01211.01211.009
observation->Enspage-2006                 12.74412.74112.71212.74112.74112.744
observation->Enspage-2007                 11.06011.06111.05211.06111.06111.060
observation->Enspage-2008                 12.10012.10412.10112.10412.10412.100
observation->Enspage-2009                 10.05510.05610.05510.05610.05610.055
observation->Enspage-2010                 17.57617.58117.62717.58117.58117.576
observation->Enspage-2011                 13.40213.39913.41213.39913.39913.402
observation->Enspage-2012                 13.14813.14913.13413.14913.14913.148
observation->Enspage-2013                 16.69116.69416.65916.69416.69416.691
observation->Enspage-2014                 13.86813.86913.87013.86913.86913.868
observation->Enspage-2015                 12.76412.76012.77612.76012.76012.764
observation->Enspage-2016                 19.22219.22419.20319.22419.22419.222
observation->Enspage-2017                 15.55215.56215.51115.56215.56215.552
observation->WnspOLF-1992                 10.12410.12210.12910.12210.12210.124
observation->WnspOLF-1993                 9.0479.0458.9919.0459.0459.047
observation->WnspOLF-1994                 15.93115.93015.89315.93015.93015.931
observation->WnspOLF-1996                 9.7849.7889.7529.7889.7889.784
observation->WnspOLF-1999                 11.55311.55311.59311.55411.55311.553
observation->WnspOLF-2000                 13.81713.81213.89213.81213.81213.817
observation->Wnspage-2001                 8.0368.0327.9948.0328.0328.036
observation->Wnspage-2002                 8.7148.7158.7018.7158.7158.714
observation->Wnspage-2003                 10.47310.47310.50710.47310.47310.473
observation->Wnspage-2004                 7.1547.1557.1407.1557.1557.154
observation->Wnspage-2006                 9.5909.5899.5519.5899.5899.590
observation->Wnspage-2007                 10.86010.86110.82210.86110.86110.860
observation->Wnspage-2008                 12.42212.42512.43512.42512.42512.422
observation->Wnspage-2009                 14.71914.71614.79014.71614.71614.719
observation->Wnspage-2010                 11.40511.40511.37011.40511.40511.405
observation->Wnspage-2011                 11.26511.26711.27511.26711.26711.265
observation->Wnspage-2012                 10.45310.45410.46510.45410.45410.453
observation->Wnspage-2013                 13.63413.63313.63413.63313.63313.634
observation->Wnspage-2014                 12.60912.60812.59712.60812.60812.609
observation->Wnspage-2016                 15.68315.68715.72215.68715.68715.683
observation->CRsumage-1992                 15.57015.58415.56615.58415.58415.570
observation->CRsumage-1993                 10.61310.61610.62410.61610.61610.613
observation->CRsumage-1994                 18.17418.17018.13018.17018.17018.174
observation->CRsumage-1995                 14.28714.28514.24214.28514.28514.287
observation->CRsumage-1996                 11.56611.56711.60711.56711.56711.566
observation->CRsumage-1997                 15.75015.75415.68815.75415.75415.750
observation->CRsumage-1998                 14.68614.69014.69014.69014.69014.686
observation->CRsumage-1999                 12.06112.06112.08412.06112.06112.061
observation->CRsumage-2000                 14.32614.32314.28314.32314.32314.326
observation->CRsumage-2001                 19.65119.65919.52019.65919.65919.651
observation->CRsumage-2002                 10.97110.96910.96010.96910.96910.971
observation->CRsumage-2003                 11.25411.25111.21011.25111.25111.254
observation->CRsumage-2004                 9.5569.5539.5619.5539.5539.556
observation->CRsumage-2005                 12.74912.75312.65412.75312.75312.749
observation->CRsumage-2006                 10.31710.31610.31710.31610.31610.317
observation->CRsumage-2007                 13.56113.56213.50413.56213.56213.561
observation->CRsumage-2008                 11.66511.66111.63911.66111.66111.665
observation->CRsumage-2009                 12.00212.00611.96512.00612.00612.002
observation->CRsumage-2010                 11.05211.04811.05511.04811.04811.052
observation->CRsumage-2011                 12.57712.57912.57112.57912.57912.577
observation->CRsumage-2012                 17.23017.22817.14117.22817.22817.230
observation->CRsumage-2013                 11.84511.84311.82211.84311.84311.845
observation->CRsumage-2014                 11.28811.28611.27111.28611.28611.288
observation->CRsumage-2016                 9.8679.8669.8479.8669.8669.867
observation->CRsumage-2018                 11.93511.93111.89511.93111.93111.935
observation->SAsumage-1992                 12.24112.24612.26912.24612.24612.241
observation->SAsumage-1993                 12.48312.47012.50212.47012.47012.483
observation->SAsumage-1994                 10.80810.81110.81810.81110.81110.808
observation->SAsumage-2001                 12.76212.74912.83912.74912.74912.762
observation->SAsumage-2002                 12.21612.20912.25412.20912.20912.216
observation->SAsumage-2003                 13.38813.38613.38613.38613.38613.388
observation->SAsumage-2004                 11.06411.06711.00511.06711.06711.064
observation->SAsumage-2005                 15.28015.27815.38615.27815.27815.280
observation->SAsumage-2006                 14.11014.11214.07814.11214.11214.110
observation->SAsumage-2007                 10.64910.64910.63910.64910.64910.649
observation->SAsumage-2008                 13.76613.77113.83613.77113.77113.766
observation->SAsumage-2009                 15.52015.51315.53015.51315.51315.520
observation->SAsumage-2010                 15.08215.08115.13015.08115.08115.082
observation->SAsumage-2012                 13.43913.43713.47913.43713.43713.439
observation->SAsumage-2013                 11.95611.95411.98211.95411.95411.956
observation->SAsumage-2015                 13.52413.52213.54013.52213.52213.524
observation->SAsumage-2017                 12.48712.49212.47812.49212.49212.487
observation->SAautage-1992                 11.81911.80911.85711.80911.80911.819
observation->SAautage-1996                 10.45010.44910.42610.44910.44910.450
observation->SAautage-1998                 12.81012.80912.83512.80912.80912.810
observation->pspawn_1993                 -6.016-6.017-6.015-6.017-6.017-6.016
observation->pspawn-1992                 -4.852-4.848-4.829-4.848-4.848-4.852
observation->pspawn-1998                 -8.291-8.284-8.300-8.284-8.284-8.291
prior->CSacousq->catchability[CSacous].q               -0.847-0.849-0.839-0.849-0.849-0.847
prior->WCacousq->catchability[WCacous].q               -1.163-1.162-1.202-1.162-1.162-1.163
prior->CRsumq->catchability[CRsum].q               -2.153-2.155-2.207-2.155-2.155-2.153
prior->SAsumq->catchability[SAsum].q               -1.760-1.762-1.557-1.762-1.762-1.760
prior->SAautq->catchability[SAaut].q               -1.911-1.913-1.762-1.913-1.913-1.911
prior->CR_process_error->observation[CRsumbio].process_error   0.0000.0000.0000.0000.0000.000
prior->SA_process_error->observation[SAsumbio].process_error   0.0000.0000.0000.0000.0000.000
prior->B0_E_with_total_log_b0_prior->process[recruit_E].b0              0.0000.0000.0000.0000.0000.000
prior->B0_W_with_proportion_prior->process[recruit_W].b0              12.54512.54412.55512.54412.54412.545
prior->M_all->process[Instant_mortality].m{west.sa}             -1.083-1.085-1.217-1.085-1.085-1.083
prior->Wspsl_shift_param->time_varying[shifted_mu].a               0.4220.4200.3780.4200.4200.422
prior->YCS_E->process[recruit_E].ycs_values{1975}-0.624-0.624-0.623-0.624-0.624-0.624
prior->YCS_E->process[recruit_E].ycs_values{1976}-0.314-0.316-0.319-0.316-0.316-0.314
prior->YCS_E->process[recruit_E].ycs_values{1977}-0.225-0.227-0.238-0.227-0.227-0.225
prior->YCS_E->process[recruit_E].ycs_values{1978}0.6560.6540.6150.6540.6540.656
prior->YCS_E->process[recruit_E].ycs_values{1979}0.0860.0840.0620.0840.0840.086
prior->YCS_E->process[recruit_E].ycs_values{1980}-0.630-0.630-0.627-0.630-0.630-0.630
prior->YCS_E->process[recruit_E].ycs_values{1981}-0.619-0.619-0.623-0.619-0.619-0.619
prior->YCS_E->process[recruit_E].ycs_values{1982}-0.551-0.550-0.542-0.550-0.550-0.551
prior->YCS_E->process[recruit_E].ycs_values{1983}-0.627-0.627-0.630-0.627-0.627-0.627
prior->YCS_E->process[recruit_E].ycs_values{1984}-0.591-0.591-0.596-0.591-0.591-0.591
prior->YCS_E->process[recruit_E].ycs_values{1985}-0.641-0.641-0.642-0.641-0.641-0.641
prior->YCS_E->process[recruit_E].ycs_values{1986}-0.408-0.408-0.400-0.408-0.408-0.408
prior->YCS_E->process[recruit_E].ycs_values{1987}0.1420.1420.1350.1420.1420.142
prior->YCS_E->process[recruit_E].ycs_values{1988}0.0250.0260.0190.0260.0260.025
prior->YCS_E->process[recruit_E].ycs_values{1989}-0.631-0.631-0.632-0.631-0.631-0.631
prior->YCS_E->process[recruit_E].ycs_values{1990}0.1590.1590.1580.1590.1590.159
prior->YCS_E->process[recruit_E].ycs_values{1991}-0.643-0.643-0.643-0.643-0.643-0.643
prior->YCS_E->process[recruit_E].ycs_values{1992}-0.643-0.643-0.643-0.643-0.643-0.643
prior->YCS_E->process[recruit_E].ycs_values{1993}-0.643-0.643-0.643-0.643-0.643-0.643
prior->YCS_E->process[recruit_E].ycs_values{1994}-0.531-0.530-0.525-0.530-0.530-0.531
prior->YCS_E->process[recruit_E].ycs_values{1995}-0.534-0.534-0.538-0.534-0.534-0.534
prior->YCS_E->process[recruit_E].ycs_values{1996}-0.573-0.574-0.577-0.574-0.574-0.573
prior->YCS_E->process[recruit_E].ycs_values{1997}-0.623-0.623-0.622-0.623-0.623-0.623
prior->YCS_E->process[recruit_E].ycs_values{1998}-0.511-0.512-0.513-0.512-0.512-0.511
prior->YCS_E->process[recruit_E].ycs_values{1999}2.0142.0142.0142.0142.0142.014
prior->YCS_E->process[recruit_E].ycs_values{2000}-0.616-0.616-0.614-0.616-0.616-0.616
prior->YCS_E->process[recruit_E].ycs_values{2001}-0.065-0.065-0.071-0.065-0.065-0.065
prior->YCS_E->process[recruit_E].ycs_values{2002}-0.542-0.543-0.545-0.543-0.543-0.542
prior->YCS_E->process[recruit_E].ycs_values{2003}-0.586-0.586-0.585-0.586-0.586-0.586
prior->YCS_E->process[recruit_E].ycs_values{2004}-0.568-0.568-0.568-0.568-0.568-0.568
prior->YCS_E->process[recruit_E].ycs_values{2005}-0.579-0.579-0.578-0.579-0.579-0.579
prior->YCS_E->process[recruit_E].ycs_values{2006}-0.483-0.483-0.482-0.483-0.483-0.483
prior->YCS_E->process[recruit_E].ycs_values{2007}-0.567-0.567-0.566-0.567-0.567-0.567
prior->YCS_E->process[recruit_E].ycs_values{2008}-0.643-0.643-0.643-0.643-0.643-0.643
prior->YCS_E->process[recruit_E].ycs_values{2009}-0.561-0.561-0.562-0.561-0.561-0.561
prior->YCS_E->process[recruit_E].ycs_values{2010}-0.031-0.031-0.046-0.031-0.031-0.031
prior->YCS_E->process[recruit_E].ycs_values{2011}-0.408-0.408-0.408-0.408-0.408-0.408
prior->YCS_E->process[recruit_E].ycs_values{2012}-0.019-0.019-0.023-0.019-0.019-0.019
prior->YCS_E->process[recruit_E].ycs_values{2013}-0.238-0.238-0.242-0.238-0.238-0.238
prior->YCS_E->process[recruit_E].ycs_values{2014}-0.317-0.317-0.314-0.317-0.317-0.317
prior->YCS_E->process[recruit_E].ycs_values{2015}-0.284-0.279-0.302-0.279-0.279-0.284
prior->YCS_E->process[recruit_E].ycs_values{2016}0.3640.3630.3730.3630.3630.364
prior->YCS_W->process[recruit_W].ycs_values{1975}-0.616-0.616-0.617-0.616-0.616-0.616
prior->YCS_W->process[recruit_W].ycs_values{1976}-0.195-0.195-0.184-0.195-0.195-0.195
prior->YCS_W->process[recruit_W].ycs_values{1977}-0.376-0.376-0.396-0.376-0.376-0.376
prior->YCS_W->process[recruit_W].ycs_values{1978}-0.397-0.396-0.417-0.396-0.396-0.397
prior->YCS_W->process[recruit_W].ycs_values{1979}-0.607-0.607-0.612-0.607-0.607-0.607
prior->YCS_W->process[recruit_W].ycs_values{1980}-0.568-0.567-0.573-0.567-0.567-0.568
prior->YCS_W->process[recruit_W].ycs_values{1981}-0.502-0.502-0.508-0.502-0.502-0.502
prior->YCS_W->process[recruit_W].ycs_values{1982}-0.606-0.606-0.607-0.606-0.606-0.606
prior->YCS_W->process[recruit_W].ycs_values{1983}-0.108-0.109-0.111-0.109-0.109-0.108
prior->YCS_W->process[recruit_W].ycs_values{1984}-0.161-0.163-0.163-0.163-0.163-0.161
prior->YCS_W->process[recruit_W].ycs_values{1985}-0.619-0.619-0.620-0.619-0.619-0.619
prior->YCS_W->process[recruit_W].ycs_values{1986}-0.097-0.098-0.098-0.098-0.098-0.097
prior->YCS_W->process[recruit_W].ycs_values{1987}0.7890.7880.7820.7880.7880.789
prior->YCS_W->process[recruit_W].ycs_values{1988}-0.251-0.252-0.255-0.252-0.252-0.251
prior->YCS_W->process[recruit_W].ycs_values{1989}1.2071.2071.2111.2071.2071.207
prior->YCS_W->process[recruit_W].ycs_values{1990}1.1981.1971.2041.1971.1971.198
prior->YCS_W->process[recruit_W].ycs_values{1991}-0.081-0.081-0.080-0.081-0.081-0.081
prior->YCS_W->process[recruit_W].ycs_values{1992}0.4800.4800.4800.4800.4800.480
prior->YCS_W->process[recruit_W].ycs_values{1993}-0.643-0.643-0.643-0.643-0.643-0.643
prior->YCS_W->process[recruit_W].ycs_values{1994}0.0140.0140.0200.0140.0140.014
prior->YCS_W->process[recruit_W].ycs_values{1995}-0.409-0.410-0.413-0.410-0.410-0.409
prior->YCS_W->process[recruit_W].ycs_values{1996}-0.301-0.301-0.306-0.301-0.301-0.301
prior->YCS_W->process[recruit_W].ycs_values{1997}-0.467-0.466-0.471-0.466-0.466-0.467
prior->YCS_W->process[recruit_W].ycs_values{1998}-0.323-0.322-0.334-0.322-0.322-0.323
prior->YCS_W->process[recruit_W].ycs_values{1999}2.0142.0142.0142.0142.0142.014
prior->YCS_W->process[recruit_W].ycs_values{2000}-0.643-0.643-0.643-0.643-0.643-0.643
prior->YCS_W->process[recruit_W].ycs_values{2001}0.8550.8560.8460.8560.8560.855
prior->YCS_W->process[recruit_W].ycs_values{2002}-0.222-0.222-0.222-0.222-0.222-0.222
prior->YCS_W->process[recruit_W].ycs_values{2003}-0.612-0.612-0.612-0.612-0.612-0.612
prior->YCS_W->process[recruit_W].ycs_values{2004}-0.570-0.570-0.571-0.570-0.570-0.570
prior->YCS_W->process[recruit_W].ycs_values{2005}-0.574-0.574-0.576-0.574-0.574-0.574
prior->YCS_W->process[recruit_W].ycs_values{2006}-0.280-0.280-0.279-0.280-0.280-0.280
prior->YCS_W->process[recruit_W].ycs_values{2007}-0.438-0.438-0.436-0.438-0.438-0.438
prior->YCS_W->process[recruit_W].ycs_values{2008}-0.464-0.464-0.463-0.464-0.464-0.464
prior->YCS_W->process[recruit_W].ycs_values{2009}-0.514-0.514-0.508-0.514-0.514-0.514
prior->YCS_W->process[recruit_W].ycs_values{2010}0.0700.0710.0780.0710.0710.070
prior->YCS_W->process[recruit_W].ycs_values{2011}0.0300.0310.0430.0310.0310.030
prior->YCS_W->process[recruit_W].ycs_values{2012}-0.049-0.050-0.059-0.050-0.050-0.049
prior->YCS_W->process[recruit_W].ycs_values{2013}-0.488-0.489-0.490-0.489-0.489-0.488
prior->YCS_W->process[recruit_W].ycs_values{2014}1.7271.7301.7421.7301.7301.727
prior->YCS_W->process[recruit_W].ycs_values{2015}1.6841.6811.6871.6811.6811.684
prior->YCS_W->process[recruit_W].ycs_values{2016}0.3640.3630.3730.3630.3630.364
prior->sel_Whome->selectivity[sel_Whome].v{1}            0.0000.0000.0000.0000.0000.000
prior->sel_Whome->selectivity[sel_Whome].v{2}            0.0000.0000.0000.0000.0000.000
prior->sel_Whome->selectivity[sel_Whome].v{3}            0.0000.0000.0000.0000.0000.000
prior->sel_Whome->selectivity[sel_Whome].v{4}            0.0000.0000.0000.0000.0000.000
prior->sel_Whome->selectivity[sel_Whome].v{5}            0.0000.0000.0000.0000.0000.000
prior->sel_Whome->selectivity[sel_Whome].v{6}            0.0000.0000.0000.0000.0000.000
prior->sel_Whome->selectivity[sel_Whome].v{7}            0.0000.0000.0000.0000.0000.000
prior->sel_Espmg_male->selectivity[sel_Espmg].v{1}            0.0000.0000.0000.0000.0000.000
prior->sel_Espmg_male->selectivity[sel_Espmg].v{2}            0.0000.0000.0000.0000.0000.000
prior->sel_Espmg_male->selectivity[sel_Espmg].v{3}            0.0000.0000.0000.0000.0000.000
prior->sel_Espmg_male->selectivity[sel_Espmg].v{4}            0.0000.0000.0000.0000.0000.000
prior->sel_Espmg_male->selectivity[sel_Espmg].v{5}            0.0000.0000.0000.0000.0000.000
prior->sel_Espmg_male->selectivity[sel_Espmg].v{6}            0.0000.0000.0000.0000.0000.000
prior->sel_Espmg_male->selectivity[sel_Espmg].v{7}            0.0000.0000.0000.0000.0000.000
prior->sel_Espmg_male->selectivity[sel_Espmg].v{8}            0.0000.0000.0000.0000.0000.000
prior->sel_Wspmg_male->selectivity[sel_Wspmg].v{1}            0.0000.0000.0000.0000.0000.000
prior->sel_Wspmg_male->selectivity[sel_Wspmg].v{2}            0.0000.0000.0000.0000.0000.000
prior->sel_Wspmg_male->selectivity[sel_Wspmg].v{3}            0.0000.0000.0000.0000.0000.000
prior->sel_Wspmg_male->selectivity[sel_Wspmg].v{4}            0.0000.0000.0000.0000.0000.000
prior->sel_Wspmg_male->selectivity[sel_Wspmg].v{5}            0.0000.0000.0000.0000.0000.000
prior->sel_Wspmg_male->selectivity[sel_Wspmg].v{6}            0.0000.0000.0000.0000.0000.000
prior->sel_Wspmg_male->selectivity[sel_Wspmg].v{7}            0.0000.0000.0000.0000.0000.000
prior->sel_Wspmg_male->selectivity[sel_Wspmg].v{8}            0.0000.0000.0000.0000.0000.000
prior->Enspsl_mu->selectivity[Enspsl].mu              0.0000.0000.0000.0000.0000.000
prior->Enspsl_s_l->selectivity[Enspsl].sigma_l         0.0000.0000.0000.0000.0000.000
prior->Enspsl_s_r->selectivity[Enspsl].sigma_r         0.0000.0000.0000.0000.0000.000
prior->Wnspsl_mu->selectivity[Wnspsl].mu              0.0000.0000.0000.0000.0000.000
prior->Wnspsl_s_l->selectivity[Wnspsl].sigma_l         0.0000.0000.0000.0000.0000.000
prior->Wnspsl_s_r->selectivity[Wnspsl].sigma_r         0.0000.0000.0000.0000.0000.000
prior->Espsl_mu->selectivity[Espsl].mu              0.0000.0000.0000.0000.0000.000
prior->Espsl_s_l->selectivity[Espsl].sigma_l         0.0000.0000.0000.0000.0000.000
prior->Espsl_s_r->selectivity[Espsl].sigma_r         0.0000.0000.0000.0000.0000.000
prior->CRsl_mu->selectivity[CRsl].mu              0.0000.0000.0000.0000.0000.000
prior->CRsl_s_l->selectivity[CRsl].sigma_l         0.0000.0000.0000.0000.0000.000
prior->CRsl_s_r->selectivity[CRsl].sigma_r         0.0000.0000.0000.0000.0000.000
prior->SAsl_mu->selectivity[SAsl].mu              0.0000.0000.0000.0000.0000.000
prior->SAsl_s_l->selectivity[SAsl].sigma_l         0.0000.0000.0000.0000.0000.000
prior->SAsl_s_r->selectivity[SAsl].sigma_r         0.0000.0000.0000.0000.0000.000
additional_prior->YCS.eq.16              0.0000.0000.0000.0000.0000.000
additional_prior->sp.migr            0.0090.0090.0130.0090.0090.009
jacobian->log_sum_b0                 0.0000.0000.0000.0000.0000.000
total_score                 1946.8001946.8101946.5401946.8101946.8101946.800

## [1] "CASAL base model convergence information"
## [1] "fmm: have converged: t = 8.24015e-007 f = 1946.54"                                                        
## [2] "Successful convergence in optimise"                                                                       
## [3] "Minimiser achieved convergence after 599 quasi-Newton iterations using 669 objective function evaluations"
## [1] ""
## [1] "CASAL sensitivity 1 model convergence information"
## [1] "fmm: step size too small. Indicates successful convergence (though this is not the textbook ideal convergence situation)"
## [2] "Successful convergence in optimise"                                                                                      
## [3] "Minimiser achieved convergence after 665 quasi-Newton iterations using 765 objective function evaluations"               
## [1] ""
## [1] "Casal2 betadiff_casal_flags_on model convergence information"
## [1] "fmm: have converged: t = 6.31891e-07 f = 1946.8"
## [2] "Successful convergence in optimise"             
## [1] ""
## [1] "Casal2 betadiff_casal_flags_off model convergence information"
## [1] "fmm: have converged: t = 4.14773e-07 f = 1946.81"
## [2] "Successful convergence in optimise"              
## [1] ""
## [1] "Casal2 betadiff_casal_flags_on_low_tol model convergence information"
## [1] "fmm: have converged: t = 6.64963e-10 f = 1946.54"
## [2] "Successful convergence in optimise"              
## [1] ""
## [1] "Casal2 adolc_casal_flags_on model convergence information"
## character(0)
## [1] ""
## [1] "Casal2 adolc_casal_flags_off model convergence information"
## character(0)
## [1] ""
## [1] "Casal2 adolc_casal_flags_on_low_tol model convergence information"
## [1] "[IMPORTANT] FMM: Small step size convergence"
## [2] "FMM: Small step size convergence"
## [1] "CASAL model warnings"
##                  parameter estimate lower.bound upper.bound
## 1   recruitment[E].YCS[25]     0.06        0.06         8.6
## 2   recruitment[W].YCS[25]     0.06        0.06         8.6
## 3 selectivity[CRsl].all[1]        1           1          17
## 4 selectivity[CRsl].all[2]        1           1          50
## [1] ""
## [1] "CASAL sensitivity 1 model warnings"
##                  parameter estimate lower.bound upper.bound
## 1   recruitment[E].YCS[25]     0.06        0.06         8.6
## 2   recruitment[W].YCS[25]     0.06        0.06         8.6
## 3 selectivity[CRsl].all[1]        1           1          17
## 4 selectivity[CRsl].all[2]        1           1          50
## [1] ""
## [1] "Casal2 betadiff_casal_flags_on model warnings"
##  [1] "Estimates were removed because of matching lower and upper bounds. Originally had 134 estimates, now have 133"
##  [2] "estimated parameter 'process[recruit_E].b0' was within 0.001 of upper bound 16.2"                             
##  [3] "estimated parameter 'process[recruit_W].b0' was within 0.001 of upper bound 0.59"                             
##  [4] "estimated parameter 'process[recruit_E].ycs_values{1999}' was within 0.001 of lower bound 0.06"               
##  [5] "estimated parameter 'process[recruit_W].ycs_values{1999}' was within 0.001 of lower bound 0.06"               
##  [6] "estimated parameter 'selectivity[CRsl].mu' was within 0.001 of lower bound 1"                                 
##  [7] "estimated parameter 'selectivity[SAsl].sigma_r' was within 0.001 of lower bound 1"                            
##  [8] "estimated parameter 'process[recruit_E].b0' was within 0.001 of upper bound 16.2"                             
##  [9] "estimated parameter 'process[recruit_W].b0' was within 0.001 of upper bound 0.59"                             
## [10] "estimated parameter 'process[recruit_E].ycs_values{1999}' was within 0.001 of lower bound 0.06"               
## [11] "estimated parameter 'process[recruit_W].ycs_values{1999}' was within 0.001 of lower bound 0.06"               
## [12] "estimated parameter 'selectivity[CRsl].mu' was within 0.001 of lower bound 1"                                 
## [13] "estimated parameter 'selectivity[SAsl].sigma_r' was within 0.001 of lower bound 1"                            
## [1] ""
## [1] "Casal2 betadiff_casal_flags_off model warnings"
##  [1] "Estimates were removed because of matching lower and upper bounds. Originally had 134 estimates, now have 133"
##  [2] "estimated parameter 'process[recruit_E].b0' was within 0.001 of upper bound 16.2"                             
##  [3] "estimated parameter 'process[recruit_W].b0' was within 0.001 of upper bound 0.59"                             
##  [4] "estimated parameter 'process[recruit_E].ycs_values{1999}' was within 0.001 of lower bound 0.06"               
##  [5] "estimated parameter 'process[recruit_W].ycs_values{1999}' was within 0.001 of lower bound 0.06"               
##  [6] "estimated parameter 'selectivity[sel_Whome].v{7}' was within 0.001 of upper bound 1"                          
##  [7] "estimated parameter 'selectivity[Enspsl].sigma_r' was within 0.001 of upper bound 50"                         
##  [8] "estimated parameter 'selectivity[Wnspsl].sigma_r' was within 0.001 of upper bound 50"                         
##  [9] "estimated parameter 'selectivity[CRsl].mu' was within 0.001 of lower bound 1"                                 
## [10] "estimated parameter 'selectivity[SAsl].sigma_r' was within 0.001 of lower bound 1"                            
## [11] "estimated parameter 'process[recruit_E].b0' was within 0.001 of upper bound 16.2"                             
## [12] "estimated parameter 'process[recruit_W].b0' was within 0.001 of upper bound 0.59"                             
## [13] "estimated parameter 'process[recruit_E].ycs_values{1999}' was within 0.001 of lower bound 0.06"               
## [14] "estimated parameter 'process[recruit_W].ycs_values{1999}' was within 0.001 of lower bound 0.06"               
## [15] "estimated parameter 'selectivity[sel_Whome].v{7}' was within 0.001 of upper bound 1"                          
## [16] "estimated parameter 'selectivity[Enspsl].sigma_r' was within 0.001 of upper bound 50"                         
## [17] "estimated parameter 'selectivity[Wnspsl].sigma_r' was within 0.001 of upper bound 50"                         
## [18] "estimated parameter 'selectivity[CRsl].mu' was within 0.001 of lower bound 1"                                 
## [19] "estimated parameter 'selectivity[SAsl].sigma_r' was within 0.001 of lower bound 1"                            
## [1] ""
## [1] "Casal2 betadiff_casal_flags_on_low_tol model warnings"
##  [1] "Estimates were removed because of matching lower and upper bounds. Originally had 134 estimates, now have 133"
##  [2] "estimated parameter 'process[recruit_E].b0' was within 0.001 of upper bound 16.2"                             
##  [3] "estimated parameter 'process[recruit_W].b0' was within 0.001 of upper bound 0.59"                             
##  [4] "estimated parameter 'process[recruit_E].ycs_values{1999}' was within 0.001 of lower bound 0.06"               
##  [5] "estimated parameter 'process[recruit_W].ycs_values{1999}' was within 0.001 of lower bound 0.06"               
##  [6] "estimated parameter 'selectivity[CRsl].mu' was within 0.001 of lower bound 1"                                 
##  [7] "estimated parameter 'process[recruit_E].b0' was within 0.001 of upper bound 16.2"                             
##  [8] "estimated parameter 'process[recruit_W].b0' was within 0.001 of upper bound 0.59"                             
##  [9] "estimated parameter 'process[recruit_E].ycs_values{1999}' was within 0.001 of lower bound 0.06"               
## [10] "estimated parameter 'process[recruit_W].ycs_values{1999}' was within 0.001 of lower bound 0.06"               
## [11] "estimated parameter 'selectivity[CRsl].mu' was within 0.001 of lower bound 1"                                 
## [1] ""
## [1] "Casal2 adolc_casal_flags_on model warnings"
##  [1] "Estimates were removed because of matching lower and upper bounds. Originally had 134 estimates, now have 133"
##  [2] "estimated parameter 'process[recruit_E].b0' was within 0.001 of upper bound 16.2"                             
##  [3] "estimated parameter 'process[recruit_W].b0' was within 0.001 of upper bound 0.59"                             
##  [4] "estimated parameter 'process[recruit_E].ycs_values{1999}' was within 0.001 of lower bound 0.06"               
##  [5] "estimated parameter 'process[recruit_W].ycs_values{1999}' was within 0.001 of lower bound 0.06"               
##  [6] "estimated parameter 'selectivity[sel_Whome].v{7}' was within 0.001 of upper bound 1"                          
##  [7] "estimated parameter 'selectivity[Enspsl].sigma_r' was within 0.001 of upper bound 50"                         
##  [8] "estimated parameter 'selectivity[Wnspsl].sigma_r' was within 0.001 of upper bound 50"                         
##  [9] "estimated parameter 'selectivity[CRsl].mu' was within 0.001 of lower bound 1"                                 
## [10] "estimated parameter 'selectivity[SAsl].sigma_r' was within 0.001 of lower bound 1"                            
## [11] "estimated parameter 'process[recruit_E].b0' was within 0.001 of upper bound 16.2"                             
## [12] "estimated parameter 'process[recruit_W].b0' was within 0.001 of upper bound 0.59"                             
## [13] "estimated parameter 'process[recruit_E].ycs_values{1999}' was within 0.001 of lower bound 0.06"               
## [14] "estimated parameter 'process[recruit_W].ycs_values{1999}' was within 0.001 of lower bound 0.06"               
## [15] "estimated parameter 'selectivity[sel_Whome].v{7}' was within 0.001 of upper bound 1"                          
## [16] "estimated parameter 'selectivity[Enspsl].sigma_r' was within 0.001 of upper bound 50"                         
## [17] "estimated parameter 'selectivity[Wnspsl].sigma_r' was within 0.001 of upper bound 50"                         
## [18] "estimated parameter 'selectivity[CRsl].mu' was within 0.001 of lower bound 1"                                 
## [19] "estimated parameter 'selectivity[SAsl].sigma_r' was within 0.001 of lower bound 1"                            
## [1] ""
## [1] "Casal2 adolc_casal_flags_off model warnings"
##  [1] "Estimates were removed because of matching lower and upper bounds. Originally had 134 estimates, now have 133"
##  [2] "estimated parameter 'process[recruit_E].b0' was within 0.001 of upper bound 16.2"                             
##  [3] "estimated parameter 'process[recruit_W].b0' was within 0.001 of upper bound 0.59"                             
##  [4] "estimated parameter 'process[recruit_E].ycs_values{1999}' was within 0.001 of lower bound 0.06"               
##  [5] "estimated parameter 'process[recruit_W].ycs_values{1999}' was within 0.001 of lower bound 0.06"               
##  [6] "estimated parameter 'selectivity[sel_Whome].v{7}' was within 0.001 of upper bound 1"                          
##  [7] "estimated parameter 'selectivity[Enspsl].sigma_r' was within 0.001 of upper bound 50"                         
##  [8] "estimated parameter 'selectivity[Wnspsl].sigma_r' was within 0.001 of upper bound 50"                         
##  [9] "estimated parameter 'selectivity[CRsl].mu' was within 0.001 of lower bound 1"                                 
## [10] "estimated parameter 'selectivity[SAsl].sigma_r' was within 0.001 of lower bound 1"                            
## [11] "estimated parameter 'process[recruit_E].b0' was within 0.001 of upper bound 16.2"                             
## [12] "estimated parameter 'process[recruit_W].b0' was within 0.001 of upper bound 0.59"                             
## [13] "estimated parameter 'process[recruit_E].ycs_values{1999}' was within 0.001 of lower bound 0.06"               
## [14] "estimated parameter 'process[recruit_W].ycs_values{1999}' was within 0.001 of lower bound 0.06"               
## [15] "estimated parameter 'selectivity[sel_Whome].v{7}' was within 0.001 of upper bound 1"                          
## [16] "estimated parameter 'selectivity[Enspsl].sigma_r' was within 0.001 of upper bound 50"                         
## [17] "estimated parameter 'selectivity[Wnspsl].sigma_r' was within 0.001 of upper bound 50"                         
## [18] "estimated parameter 'selectivity[CRsl].mu' was within 0.001 of lower bound 1"                                 
## [19] "estimated parameter 'selectivity[SAsl].sigma_r' was within 0.001 of lower bound 1"                            
## [1] ""
## [1] "Casal2 adolc_casal_flags_on_low_tol model warnings"
##  [1] "Estimates were removed because of matching lower and upper bounds. Originally had 134 estimates, now have 133"
##  [2] "estimated parameter 'process[recruit_E].b0' was within 0.001 of upper bound 16.2"                             
##  [3] "estimated parameter 'process[recruit_W].b0' was within 0.001 of upper bound 0.59"                             
##  [4] "estimated parameter 'process[recruit_E].ycs_values{1999}' was within 0.001 of lower bound 0.06"               
##  [5] "estimated parameter 'process[recruit_W].ycs_values{1999}' was within 0.001 of lower bound 0.06"               
##  [6] "estimated parameter 'selectivity[Enspsl].sigma_r' was within 0.001 of upper bound 50"                         
##  [7] "estimated parameter 'selectivity[Wnspsl].sigma_r' was within 0.001 of upper bound 50"                         
##  [8] "estimated parameter 'selectivity[CRsl].mu' was within 0.001 of lower bound 1"                                 
##  [9] "estimated parameter 'selectivity[SAsl].sigma_r' was within 0.001 of lower bound 1"                            
## [10] "estimated parameter 'process[recruit_E].b0' was within 0.001 of upper bound 16.2"                             
## [11] "estimated parameter 'process[recruit_W].b0' was within 0.001 of upper bound 0.59"                             
## [12] "estimated parameter 'process[recruit_E].ycs_values{1999}' was within 0.001 of lower bound 0.06"               
## [13] "estimated parameter 'process[recruit_W].ycs_values{1999}' was within 0.001 of lower bound 0.06"               
## [14] "estimated parameter 'selectivity[Enspsl].sigma_r' was within 0.001 of upper bound 50"                         
## [15] "estimated parameter 'selectivity[Wnspsl].sigma_r' was within 0.001 of upper bound 50"                         
## [16] "estimated parameter 'selectivity[CRsl].mu' was within 0.001 of lower bound 1"                                 
## [17] "estimated parameter 'selectivity[SAsl].sigma_r' was within 0.001 of lower bound 1"

1.5 Matching of outputs

Time series comparisons with CASAL base model results

## [1] "Catch time series base model comparison for run betadiff_casal_flags_on"
## [1] "Actual catches for Ensp1 match: yes"
## [1] "Actual catches for Wnsp1 match: yes"
## [1] "Actual catches for Ensp2 match: yes"
## [1] "Actual catches for Wnsp2 match: yes"
## [1] "Actual catches for Esp match: yes"
## [1] "Actual catches for Wsp match: yes"
## [1] ""
## [1] "Catch time series base model comparison for run betadiff_casal_flags_off"
## [1] "Actual catches for Ensp1 match: yes"
## [1] "Actual catches for Wnsp1 match: yes"
## [1] "Actual catches for Ensp2 match: yes"
## [1] "Actual catches for Wnsp2 match: yes"
## [1] "Actual catches for Esp match: yes"
## [1] "Actual catches for Wsp match: yes"
## [1] ""
## [1] "Catch time series base model comparison for run betadiff_casal_flags_on_low_tol"
## [1] "Actual catches for Ensp1 match: yes"
## [1] "Actual catches for Wnsp1 match: yes"
## [1] "Actual catches for Ensp2 match: yes"
## [1] "Actual catches for Wnsp2 match: yes"
## [1] "Actual catches for Esp match: yes"
## [1] "Actual catches for Wsp match: yes"
## [1] ""
## [1] "Catch time series base model comparison for run adolc_casal_flags_on"
## [1] "Actual catches for Ensp1 match: yes"
## [1] "Actual catches for Wnsp1 match: yes"
## [1] "Actual catches for Ensp2 match: yes"
## [1] "Actual catches for Wnsp2 match: yes"
## [1] "Actual catches for Esp match: yes"
## [1] "Actual catches for Wsp match: yes"
## [1] ""
## [1] "Catch time series base model comparison for run adolc_casal_flags_off"
## [1] "Actual catches for Ensp1 match: yes"
## [1] "Actual catches for Wnsp1 match: yes"
## [1] "Actual catches for Ensp2 match: yes"
## [1] "Actual catches for Wnsp2 match: yes"
## [1] "Actual catches for Esp match: yes"
## [1] "Actual catches for Wsp match: yes"
## [1] ""
## [1] "Catch time series base model comparison for run adolc_casal_flags_on_low_tol"
## [1] "Actual catches for Ensp1 match: yes"
## [1] "Actual catches for Wnsp1 match: yes"
## [1] "Actual catches for Ensp2 match: yes"
## [1] "Actual catches for Wnsp2 match: yes"
## [1] "Actual catches for Esp match: yes"
## [1] "Actual catches for Wsp match: yes"
## [1] ""

Derived quantities

SB0, SBcurrent, MSY, F_MSY, others…

1.6 Plots

Comparison plots

## [1] "CASAL base parameter correlation range (excluding 1.0): -0.9914 0.9969"

## [1] "CASAL sensitivity 1 parameter correlation range (excluding 1.0): -0.9324 0.9968"
## Warning in min(x): no non-missing arguments to min; returning Inf
## Warning in max(x): no non-missing arguments to max; returning -Inf

## [1] "Casal2 betadiff_casal_flags_on parameter correlation range (excluding 1.0): Inf -Inf"
## Warning in min(x): no non-missing arguments to min; returning Inf

## Warning in min(x): no non-missing arguments to max; returning -Inf

## [1] "Casal2 betadiff_casal_flags_off parameter correlation range (excluding 1.0): Inf -Inf"

## [1] "Casal2 betadiff_casal_flags_on_low_tol parameter correlation range (excluding 1.0): -0.932326 0.99673"

## [1] "Casal2 adolc_casal_flags_on parameter correlation range (excluding 1.0): -0.885609 0.996665"

## [1] "Casal2 adolc_casal_flags_off parameter correlation range (excluding 1.0): -0.888378 0.996666"

## [1] "Casal2 adolc_casal_flags_on_low_tol parameter correlation range (excluding 1.0): -0.993528 0.99951"